Plane waves as tractor beams 
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It is shown that in a large class of systems plane waves can act as tractor beams: i.e., an incident 
plane wave can exert a pulling force on the scatterer. The underlying physical mechanism for the 
pulling force is due to the sufficiently strong scattering of the incoming wave into another mode 
having a larger wave number, in which case excess momentum is created behind the scatterer. Such 
a tractor beam or negative radiation pressure effect arises naturally in systems where the coupling 
between the scattering channels is due to Aharonov-Bohm (AB) gauge potentials. It is demonstrated 
£f} . that this effect is also present if the AB potential is an induced, ("artificial") gauge potential such 

as the one found in J. March-Russell, J. Preskill, F. Wilczek, Phys. Rev. Lett. 58 2567 (1992). 
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a. 

Introduction It seems to be a hitherto unremarked phenomenon, that in multi-channel scatterings the force acting 
on the scatterer can be counter-intuitive. We show that for a large class of scattering systems with at least two channels 
with different dispersion relations, an incoming plane wave can exert a pulling force on the scatterer (tractor beam 
effect or negative radiation pressure, NRP). Models where negative radiation pressure appears have many physical 
applications, ranging from various condensed matter systems with vortices to cosmic strings. 

We shall consider in more detail systems with vorticity, where the coupling between the two channels is due to an 
^7* , Aharonov-Bohm (AB) potential. Typical cases of such systems are vortices ,where the AB potential is induced by 
Qh the "internal frame dragging" due to the breaking of a global U(l) symmetry For a recent review on systems with 
similar "artificial gaugepotential" we refer to the review 0|. Another noteworthy example is the scattering of sound 
' '■ on a superfluid vortex 0,01 or the scattering of various fields on cosmic strings [H, [H, [fj. Our prototype example will 
be the model studied in Ref. [lj], which describes the scattering of scalar fields on a global vortex or that of neutral 
particles on a magnetic vortex. Unlike Ref. [lj], we consider the case when both channels are open. We show, that for 
a large range of the parameters of the model, the incident plane wave acts as a tractor beam. 

Negative radiation pressure on localized solutions (kinks) due to an incoming plane wave has already been observed 
in one dimensional classical field theories 0, S|- I n that case NRP is induced in systems with a small reflection 
coefficient and sufficiently strong conversion to higher harmonics by the nonlinearities 0|- More recently, there has 
£f) \ been an upsurge of interest for tractor beams in optics, with light beams specifically prepared for the scatterer in a 
■ way that it is scattered dominantly in the direction behind the scatterer [9(. 

Setup In this paper, we consider stationary scattering, in one and two dimensional cases. The incoming plane 
wave propagates in the positive direction along the x axis, with the scatterer placed at the origin. We assume two 
channels u, d, with wave numbers k u , kd, related to the wave frequency ui, by the dispersion relations k 2 + m 2 = oj 2 . 

One dimensional case For the case of an incoming wave in the u channel with amplitude A u and wave number 
k u , the force acting on the scatterer follows from momentum conservation: 
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F U /A 2 U = k u {l+ \R UU \ 2 - \T UU \ 2 ) + k d {\R du \ 2 - \T du \ 2 ) , (1) 

where Rij resp. Ty denote the transmission resp. reflection coefficients for an incoming wave of type j into channel i, 
and k d is the wave number in the d channel. Eq. ([TJ) follows immediately by taking into account that a plane wave 
Aexp(ikx)/ vk, carries momentum fc|A| 2 . We note that the energy flux of such a wave is w|A| 2 . 

The conservation of energy implies that the energy flux of the outgoing waves is equal to that of the incoming wave, 
therefore 



\T ij \ 2 ) = l, (2) 



which corresponds to the unitarity of the S-matrix. By combining eqns. ([lj and ©, the first term in eq. ([T]) is positive, 
and thus for a single channel scattering Eq.([T]) reduces to the familiar radiation pressure formula F = 2/c|i?| 2 > 0. 

Already when the coupling of the channels is due to a simple step potential, the forward scattering amplitude from 
channel u to d, Td u , can become large enough to reverse the direction of the force ([1} acting on the scatterer, provided 
that kd > k u , (which is the case when, e.g., m u > md), that is NRP is induced by flavor changing scattering. 
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Physically, such a scenario can be realized by the low energy scattering of neutral particles with magnetic moment, 
in a magnetic field perpendicular to their propagation (e.g., B z < 0, the two species being spin up/down). The 
scatterer is a region with magnetic field along the x axis. 

We have illustrated that NRP on a scatterer appears for various one dimensional systems. The physical mechanism 
for the NRP in this simple setting is that if an incoming wave in the u channel is scattered sufficiently strongly to the 
d channel with > k u , excess momentum is created behind the scatterer, since the scattered wave in the d ("light") 
channel carries more momentum than the incoming one. 

Force acting on the scatterer in planar problems As we shall show, NRP occurs for multi-channel scattering in a 
number of 2 dimensional systems, in spite of waves being scattered in all directions, not just forward and backwards 
as in 1 dimension. Exploiting the conservation equations V = — VT where V is the momentum density of the field, 
we obtain the force acting on the scatterer in two dimensions: 

F = lim / Te r Rd$ , (3) 

J-n 

where T is the stress tensor, the overbar denotes time averaging, and r, $ denote the polar coordinates of the plane. 
To compute T asymptotically, we assume that in the far field region the waves can be considered freely propagating. 
Assuming n channels, for each partial wave £, the S- matrix elements form an n x n unitary matrix Si. Evaluating 
the expression ([3]) we obtain a master equation for the force F, in terms of Sf. 

F = F x + iF y = -4 J2 {^Sl +1 KS e A - A^KA} , (4) 
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where A = (A\, . . . , A n ) T , K = diag(fci, . . . , k n ) and f denotes the adjoint (transposed conjugate). Eq. (|4]) can also 
derived for potential scattering in quantum mechanics (in h — 1 units) . 

Analogously to the one dimensional case, it is not difficult to see that for a free, two component wave (it, <i), coupled 
in a disk- like region r < L, by a constant potential, NRP is observed for various masses of the particles in a large 
range of frequencies. We should like to mention that the force acting on the solenoid for the Aharonov-Bohm (AB) 
scattering , follows from our general result Eq. (j4|) . To show this we recall that the scattering phases in the partial 
wave expansion for the AB problem are given as Si — exp(2iSi), Sg — tt(£ — v)/2, with v 2 = (£ — $) 2 , where $ is the 
magnetic flux in the solenoid. Then from Eq. the well known result for the longitudinal and transversal component 



of the force [llj follows: 

F x +iF y = -A\A\ 2 k{e- 2 ^ -1). (5) 



It is a somewhat striking result that the force acting on the solenoid has also a transversal component. As it is clear 
by now this is due to the breaking of the reflection symmetry with respect to the x axis by the magnetic flux (direction 
of circulation) . 

The result ([5]) has other applications as well. For example the force exerted on a superfluid vortex by the scattering 
of sound waves in the Gross-Pitaevskii (GP) model is given by Eq. |5J with a suitable replacement of the flux 
0. The longitudinal component, F x , is the acoustic drag, and the transversal component is known as the Iordanskii 
force 0| . 

Scattering on a vortex In Ref. [l| it has been shown that in the low momentum transfer limit, the scattering cross 
section of a complex scalar ip 7 on a vortex in a GP-type model is very close to that of the AB problem. In Ref. [l| 
the problem is reduced to a single channel scattering, which is valid below the threshold defined by the larger mass 
eigenvalue m u . In the following we shall show that at higher energies when both channels are open, the scattering of 
the u mode induces NRP on the vortex. The coupling term between tp and the symmetry breaking complex scalar 
field (f>, is given by AC = g(f>ip 2 + c.c. The field equations for the field if) are 

(uj 2 - m 2 + V 2 )il) = 2g(/)*^* . (6) 

Far from the vortex <f> ~ vexp(i$), therefore the mass eigenstates of the potential are defined by the transformation 
(ip,ip*) T = y/2Up, p — (u,d) T . U is a unitary matrix given as 

V2U=l„ /a _.^ /2 \ . (7) 
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Transformation yields a heavy mode u 
written as 



2gv) and a light mode d {jn 2 d 



2gv). Eqs. © can be 



(v+, A? y 



p-K*p = 0, a 2 = 



(8) 



where A = e#/r is an induced, "artificial" gauge potential with vorticity (or winding number) ±1/2, e# is the angular 
unit vector. Note, that the exp(i#/2) terms in Eq. © induce anti-periodic boundary conditions on p. 
Partial waves are introduced taking into account the vortex number of A, as 



CM) 



(9) 



with 7 = 1/2. The radial functions satisfy 
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where rfc = rfc = (£ + l/2) 2 + l/4 and c = i(£ + 1/2). We use rescaled variables such that v = 1, m u = 2 (threshold 
at u) = 2), and for simplicity sake, present numerical data for = 1. 

We have computed the S-matrix elements by solving numerically Eqs. (JTUJ) for < I < 12 for various values of the 
frequency 2 < uj < 4 and evaluated the force acting on the (point)vortex from the master Eq. (j4|), comp. Figs. 1-3. 
(We have verified that in this frequency range partial waves with I > 12 give negligible contribution.) To understand 
the results for the coupled channel scattering better, we present some analytical results based on perturbation theory 
in the coupling, c, near the threshold of the u-channel (k u <C 1). We start at the threshold (k u = 0) when channel 
u is still closed while channel d is open. Then the 0-th order solutions are: u f ] = 0; df ] = i^J Vd (k d r), with 
£d = (£ — r)d)Tr/2. Furthermore not only the first order correction can be calculated in closed form, (d^ = 0) 

(2) 

but also the asymptotic form of the second order one, d\ . The ^-matrix is reduced to a single scattering phase Sd, 
given as: 



5d = £d + A , with A = arctan 



2 tt/4 
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What one learns from the above is that Sd at the threshold is determined not just by the expected AB phase shift 
corresponding to $ = 1/2 but that the 2nd order correction, A, is still important for not too large values of £. The 
perturbative computation of the S'-matrix above the threshold is significantly simplified by taking into account the 
above phase correction, Eq. (fTTj) . for the 0-th order d channel wave, i.e., df 1 ^ = i l e l5d J L , d (kdr), where = I — 28d/^- 
Then it is sufficient to go to just 1st order in perturbation theory to obtain the leading part of the S'-matrix close to 
the threshold. We find 



1 / e 215 " 
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with / given as 



k d 



^(E/^^ + l^/fcg) 

2r(i-z/)i>„ + i) 



(13) 



where v u — r\ Ul S u = ir(£ — v u )/2), v — (y u + Vd)/2 1 v' — (y u — Vd)/2 1 and N is a normalization factor to achieve 
exact unitarity in each order. Our perturbative results provide a satisfactory qualitative description of the numerical 
solution of Eqs. (fTU)) . comp. Figs. 1-3. We remark also that this perturbative computation induces relatively large 
errors for small values of £, caused by the behavior at r = of the phase corrected c?-channel wave. 

Force acting on the scatterer Once the S-matrix is known, our master formula Eq. (j4|) yields the force acting on the 
scatterer. Consider first the decoupled problem, when the S-matrix elements are very simple, Si — diag(e 2i £" , e 2 ^ d ), 
where t; u- d — ft{£ — f]u,d)- For the longitudinal component one finds: F x ,i ~ 8.67fci (there is no induced NRP) and 



4 



+ 



4 i fi Q 



+ : □ 



q + + 



Re S uu + 

Im S uu X 

Re S uu p □ 

Im S uu p O 



Figure 1. Real and imaginary parts of S uu for oj = 2.22 in function of £, computed numerically and perturbatively (p). 




x x 



x x 





+ 


F x.d 


X 






F x,d 






F xu p 


□ 


F x/ 


o 



2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 



Figure 2. The longitudinal force component F x as a function of the frequency uj: the numerical result is compared to the 
decoupled, (dc), resp. to the perturbative, (p), approximation, (dc) resp. (p) are depicted for 2 < uj < 2.3 resp. 2 < uj < 2.5. 
For the incoming heavy mode u, the radiation pressure becomes negative at to m 2.1557. 
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Figure 3. The numerical and the perturbative (p) result for the transversal force component F y as function of the frequency uj. 
(p) is depicted only for 2 < uj J$ 2.5). 



F y = (there is no transversal force), comp. Fig. [5] Note that the drag force due to the incoming d mode is quite 
well reproduced in the decoupled approximation, while it is completely wrong for the u-mode. 
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Figure 4. The effect of the vortex core: the longitudinal force F x acting on a point vortex and on a vortex with a linear core 
(marked with 'c'). The existence of a vortex core enhances NRP. 



Slightly above the threshold, the scattering from channel u to d (the Sdu term) can lead to NRP on the vortex, 
comp. Fig. 2. The transversal force acting on the scatterer, F y , is an analog of the IordansMi force acting on vortices 
in superfiuids In our case the total transversal force is entirely due to the coupling between the channels, comp. 

Pig. 01 

As Figs. O [3l illustrate, near threshold perturbation theory describes qualitatively the effects of the coupling between 
the two modes, losing its validity for frequencies u > 2.4. The main source of errors in our perturbative approach 
appears in the phases of the matrix elements of S u d an d Sdu- 

In the preceding analysis, the force acting on a scatterer was obtained from a stationary waveform. A condition for 
the validity of this description is that the outgoing wave packets in both modes reach the asymptotic region (where 
the wave field can be described by the 5-matrix) within a characteristic time scale, which is problematic precisely 
at the threshold. For this reason, it is important for the validity of our perturbative approach that the perturbative 
results are valid for energies away from the threshold. 

Effect of a vortex profile Up to now we have considered point vortices, when the field configuration is approximated 
everywhere by its asymptotic form 4> = v exp(-ii9). A more realistic vortex configuration, cp = f(r) exp(i$), has of course 
a non-singular core region which can be described by a profile function /, satisfying /(0) = and f(r — > oo) = v. To 
illustrate how a vortex core affects our results we consider a simple profile function f(r < R c ) — vr/R c , f(r > R c ) = v. 
In Fig. HI the effect of a linear core on the longitudinal force is depicted, for R c = 1. It can be seen that the NRP 
effect is enhanced by the profile. The reason for this enhancement can be understood in that a localized potential 
affects significantly partial waves with small values of £, and brings the £ — component of 5 U closer to that of I = 1, 
resulting in the reduction of the t — —1,0 force components. 

Generalizations A remarkably large class of scattering problems, like the perturbations of superfluid vortices [3|, 
perturbations of global cosmic strings (H, 0] or the scattering of neutral spin-1/2 particles on a magnetic vortex [lj 
can be brought to the form of a multi-component AB-type scattering, with a potential matrix 

'-(,£- "ST <"> 

Diagonalizing P with a unitary matrix U and introducing mass eigenstates with tpi — \/2Uijpj, we assume field 
equations of the form 

(V - iTKfp + K 2 p = 0, T=h . (15) 

When partial waves are introduced, the coefficients of the 1/r 2 term are given as rjf = (I — 7 — e^) 2 + \q\ 2 and 
c = q(2(£ + 7) — (ei + e2)). From this point on the perturbative calculations can be carried out in the same way as 
for Eqs. p0|, 

A number of important problems can be considered in the framework of (|15j) . They include the example considered 
in Ref. [l|, the scattering of neutral spin-1/2 particles on a magnetic vortex, when q = i/2, k\ — fc 2 = gM, M being 
the magnetization. Another noteworthy case is the scattering of perturbations on superfluid vortices in the GP model 
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with 1st order dynamics when e% — — e2 = — cj/2 + q = 1, fc 2 = 2 — y/A + w 2 , A: 2 = 2 + vTTcJ 2 , (note that the 
u channel is closed for any value of to). Vortex solutions in the GP model with 2nd order dynamics are used to model 
(straight) global cosmic strings, and scattering of various fields on them has been intensively studied [H|. Scattering 
of scalar perturbations on cosmic strings fits naturally in our framework too, e\ = e2 = 0, q = — 1, k\ = to 2 — 4, 
kd = to. We have carried out a detailed study of this problem with the result that by the scattering of the u-mode 
NRP appears Q . In this case the effect of the vortex core is crucial for the NRP, which is absent in the point vortex 
limit. We also note that for the global cosmic string, due to the symmetry of the radial equations, I — > —£, the 
transversal force is absent, F y = 0. We have been able to confirm the validity of our computations by carrying out 
time-dependent numerical simulations of the field equations 
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